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INTRODUCTION 


Two projects are summarized in this report. The first project is entitled Stiffener Crippling Ini- 
tiated by Delamination , and the second is entitled Pressure Pillowing of an Orthogonally Stiff- 
ened Cylindrical Shell. 

The write-up in the “Research Accomplished” section for the first project is the text of a con- 
ference paper that appeared in 

• Mechanics of Composite Materials— Nonlinear Effects, edited by M. W. Hyer, AMD-Vol. 
159, American Society of Mechanical Engineers, New York, NY, 1993, pp. 19-28. 

This conference paper was submitted March 1, 1993 to the editor. The papers in this ASME vol- 
ume were presented at the 1st Joint Mechanics Meeting of ASME • ASCE • SES - MEET’ N’ 93, 
Charlottesville, Virginia, June 6 -9, 1993. 

The write-up in the “Research Accomplished” section for the second project is based on the 
presentation (presenter in boldface font) 

• Johnson, Eric R., and Rastogi, Naveen, “Interacting Loads in an Orthogonally Stiffened 
Composite Cylindrical Shell,” presented in Session 25 — Work-in-Progress 1, at The 34th 
AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, 
April 19-21, 1993, Hyatt Regency, LaJolla, California 


RESEARCH ACCOMPLISHED 


I. Stiffener Crippling Initiated by Delamination 

Background and objective 

Graphite-epoxy stiffeners tested in compression (Bonanni, Johnson, and Starnes, 1991) 
showed evidence of failure initiation in postbuckling by delamination at the free edges of the 
flanges. This mode of failure initiation occurred in I- and /-section specimens with flange width to 
thickness ratios in the range from ten to twenty. The analyses conducted in the Bonanni study 
modeled the stiffeners as a branched shell using the STAGS computer code. Since shell theories 
are based on plane stress and neglect through-the-thickness stress components, failure initiated by 
interlaminar stresses at the free edges of the flanges due to bending in postbuckling could not be 
studied with the STAGS model. 

The objective of this project is to develop a computational model of the stiffener specimens 
that includes the capability to predict the interlaminar stress response at the flange free edge in 
postbuckling. 
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Finite Element Analysis 


Variable complexity modeling 

The analyses were performed using the COMET (Stewart, 1989) structural analysis software 
system. A variable complexity modeling strategy is employed which first captures the overall geo- 
metrically nonlinear response of the specimen to its failure initiation load using all shell elements 
in the model, and second changes the complexity of the model to include a sub-domain of solid 
elements at the free edge of a flange from which interlaminar stresses can be computed. The 
branched shell model of the stiffener is shown in Fig. 1, and the complex model is shown in Fig. 
2. The displacement field of the shell model is used to estimate the initial guess of displacement 
field for the complex model in Newton’s method at the failure load. Thus, the complex model only 
goes through one sequence of Newton’s iterations at the final load step, rather than using the com- 
plex model through all the intermediate load steps. This strategy keeps problem size small relative 
to a fully nonlinear three-dimensional model of the I-stiffener. A fully three-dimensional model 
would tax available computer resources and is not warranted based on the localized nature of the 
three-dimensional effects. 

Shell, transition, and solid elements 

The shell element used from the COMET library is denoted EX97, and it is an assumed natu- 
ral strain element with nine nodes and five degrees of freedom per node. The degrees of freedom 
are the three displacements of the reference surface and the two rotations of the normal to the ref- 
erence surface. (No drilling degree of freedom is included.) Element EX97 can be used for flat 
plate analysis, which is the situation for the web and flanges of the I-stiffener. The solid element is 
the standard twenty-node serendipity brick with three displacement degrees of freedom per node. 
This brick element is denoted as BR20 in the COMET library. Since the shell and brick elements 
are kinematically incompatible, transition elements are required to interface between them. The 
two transition elements TR12 and TR15 in COMET are employed (Dtivila, 1993). Element TR 12 
has twelve nodes, five of which have the five shell degrees of freedom compatible with EX97 , and 
the remaining seven nodes have the three displacement degrees of freedom compatible with 
BR20. Element TR15 has fifteen nodes, three of which have the five shell degrees of freedom and 
the remaining twelve nodes have the three solid element degrees of freedom. These elements are 
particularly useful for laminate modeling since one shell element can be connected to several 
TR12 or TR15 elements stacked through the thickness of the laminate. That is, the nodes of the 
transition element that have the five shell degrees of freedom do not have to lie on the boundary of 
the element. Element TR12 is used at two of the comers of the sub-domain since it is compatible 
on two of its adjacent faces with EX97, whereas element TR15 has one shell compatible face and 
is used along three of the edges of the sub-domain. BR20 elements fill the interior of the sub- 
domain and the interior part of the fourth edge, which coincides with the traction-free edge for the 
flange. 

The geometrically nonlinear analysis is based on the total Lagrangian formulation, and the 
fully nonlinear Green-Lagrange strain-displacement relations are employed for the solid and tran- 
sition elements. A modified Newton’s method is used to solve the nonlinear equations, and an arc- 
length strategy is used for continuation along the nonlinear equilibrium path in load-displacement 
space. 

Modeling Specimens 12 and 110 

Analyses of specimens 12 and 110 from the experimental study by Bonanni, et al. (1991) were 
conducted using the methodology discussed above. Nomenclature for the cross-sectional dimen- 
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sions is shown in Fig. 3., and measured dimensional data used in the analysis of each specimen 
arc listed in Table 1. Also, the buckling load P CT and crippling load P cc from the experiment are 
listed in this table. The specimens were fabricated from AS4/3502 graphite/epoxy unidirectional 
tape and material properties used are listed in Table 2. The wall laminate for specimen 12 was an 
eight ply layup with a stacking sequence of [±45/0/90] s , where a zero degree ply designates fibers 
parallel to the x-axis in Fig. 1. The laminate for specimen 110 is sixteen plies with [±45/0/90] 2 s 
stacking sequence. However, the fabrication procedure for the 1-specimens resulted in two diago- 
nally opposite flanges having unsymmetrical stacking sequences. For die eight-ply wall the 
unsymmetrical stacking sequence is [±45/0/902/0/±45 ]j, and for the sixteen-ply wall the unsym- 
metrical stacking sequence is [(±45/0/90)2/(90/0/±45)2 ]t- About one inch of the stiffener at each 
end was supported in an aluminum potting compound as part of the end fixture support. The 
unsupported length of specimen outside the potted ends is denoted as the gage length lg, and the 
total length of the specimen is denoted by L. Clamped end boundary conditions were imposed at 
each end (x = 0 and x = L) in the finite element models, as well as prescribing the out-of-plane 
displacement for the elements contained in the potted ends to vanish. 

Credibility of the branched shell model for each stiffener was based on comparing the finite 
element analysis to the experiment on the load-end shortening response plot. This plot for speci- 
men 110 is shown in Fig. 4. (The correlation between analysis and experiment for specimen 12 
was better than shown for 110 in Fig. 4.) To achieve the good agreement between analysis and 
experiment as shown in Fig. 4, we determined by trial and error selected material stiffnesses to 
closely match the prebuckling stiffness and buckling load of the specimens in the experiment. The 
unidirectional filler material in each web-flange junction shown in Fig. 3 was modeled as a beam- 
column and the modulus selected for this material is listed in Table 2. In addition, the fiber direc- 
tion modulus for AS4/3502 was changed from its nominal value of 18.5 Msi, and the value 
selected for each specimen is listed in Table 1. 

There were 7,175 degrees of freedom for the unrestrained branched shell model of specimen 
12, and 5,945 degrees of freedom for the unrestrained branched shell model of specimen 110. The 
complex model for specimen 12 had 168 BR20 elements and 10,325 degrees of freedom, and the 
complex model for specimen 110 had 192 BR20 elements and 9,005 degrees of freedom. The sub- 
domain of solid elements was about one inch in length (x-direction in Fig. 1) and its position 
along the flange could be varied. The width (y-direction) of the sub-domain was 0.30 inches for 
specimen 12, or 4.4/ where / is the laminate’s thickness. The width was 0.20 inches, or 2.4/, for 
specimen 110. Note that coordinate y = 0 at the free edge. The number of solid elements through 
the thickness of the flange for specimens 12 and 110 is shown in Fig. 5. When one element is used 
to span more than one ply, then effective three-dimensional stiffnesses for this sublaminate are 
determined by a homogenization procedure contained in the LAU processor of COMET. This 
homogenization procedure is similar to one presented by Sun and Li (1988). 

Delamination Index 

The intensity of the interlaminar stresses was quantified by the dimensionless index F, defined 
by 
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where x^ is the interlaminar normal stress (positive in tension), t h is the interlaminar shear stress 
tangent to the free edge, and x zy is the interlaminar shear stress normal to the free edge. The inter- 
laminar strengths Z T , Z S1 , and Z S2 used for the AS4/3502 graphite/epoxy material system are 
given in Table 2. Since this data is difficult to measure, the interlaminar tensile strength is 

assumed equal to the lamina in-plane transverse tensile strength (Z T = Y{), and the interlaminar 

shear strengths are assumed equal to the lamina in-plane shear strength (Z = Z — S). Only ten- 
sile interlaminar normal stresses are assumed to contribute to the delamination index. The form of 
Eq. (1) suggests that values of F> 1 define an interlaminar strength failure. However, it is likely 
that an elasticity solution to this effective modulus model of the laminate would have a stress sin- 
gularity at the free edge, and this implies a criterion based on stress magnitude is meaningless. 
Thus, we prefer to use Eq. (1) as an indicator of the severity of the interlaminar stress state for this 
approximate solution (finite element model). 

Equation (1) is similar to the quadratic equations used for the prediction of delamination initi- 
ation at the straight free edge in tensile/compression coupons by Brewer and Lagace (1988) and 
by Kim and Soni (1986). These delamination criteria use average values of the interlaminar 
stresses rather than point-wise values. Average values of the interlaminar stresses are obtained by 
integrating the interlaminar distributions over a critical length in the interface, and this length is 
hypothesized to be one ply thickness or a material parameter. Also, the manner in which compres- 
sive interlaminar normal stress contributes to these quadratic stress criteria are different than what 
is shown in Eq. (1). 


Results And Discussion 

The distribution of the delamination index for specimen 110 is plotted over the 0/90 interface 
in Fig. 6 at a load near failure (P/P cc = 1.03). This 0/90 interface is near the mid-plane of the lam- 
inate and is labeled number 3 in Fig. 5(b). The longitudinal position of the sub-domain of solid 
elements is 1.6 in. < x < 2.4 in. This location corresponds to the neighborhood of the inflection 
point in the plot of the out-of-plane displacement of the flange along its free edge from the shell 
solution, w(x, 0), as shown in Fig. 7. The largest magnitude of Fin Fig. 6 is 0.73 (= Fj^) and it 
occurs at the free edge where y = 0 (actually it occurs at the Gauss point closest to the free edge) 
and at x = 2.03 in. This data is also itemized in Table 3. The largest contributor to F^ is the 
interlaminar shear stress as shown in Table 3, and the interlaminar normal stress is compres- 
sive at this location. The distributions of x zx and x zz in the 0/90 interface are shown in Figs. 8 and 
9, respectively. Other longitudinal positions for the sub-domain were investigated, but the neigh- 
borhood of the inflection point resulted in the largest value of F found. For instance, the maximum 
value of F in the neighborhood of the buckle crest is .0.02 (see Table 3). 

For specimen 12 at a load near failure (F/P c c = 1 02) F m ^ = 0.4, and it occurs at x = 6.24 in., y 
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= 0, and in the 4th interface (90/0) shown in Fig. 5(a). The longitudinal position where F^^ 
occurs is again near an inflection point of w(x,0) from the shell solution. The largest contributor to 
F aax is the interlaminar shear stress x zx , and x zz is compressive. Longitudinal positions near the 
buckle crest and node point were investigated as well, and these results are listed in Table 3. 

In both specimens the severest interlaminar stress intensity is dominated by the shear stress 
component (x M ) tangent to the free edge. The axial location along the free edge of the flange 
where Fmax occurs is near the location where the transverse shear resultant Q x from the shell solu- 
tion is a relative maximum. Refer to Fig. 7. Relative maxima of Q x are near the inflection points 
of w(x,0). Thus, x w attains a large value where its corresponding shell resultant Q x is large, which 

is what one would expect. For postbuckled laminated plates having edge support all around Star- 
nes and Rouse (1981), and Davies et al. (1986), have also observed that the transverse shear 
resultant, or the shear stress component tangent to the edge, causes failure at the intersection 
nodal line and the unloaded edge. The interlaminar normal stress component (x^ is compressive 
where occurs, and its magnitude (Table 3) is small; i.e., the square of the ratio of x^ to the 

transverse compression strength of a lamina, which is about 30 ksi for AS4/3502, is small com- 
pared to the shear terms in Eq. (1). Thus, for the analyses of specimens 12 and 110 it appears 
unlikely that the compression values of x^ would contribute much to delamination initiation. 

We expect the maximum value of the index F to exceed unity if delamination initiation is pre- 
dicted. Failure initiation for specimen 110 was observed in the experiment to be by free edge 
delamination, and an F^ = 0.73 at P = 1.03 P cc from the analysis does not correlate with this 

observation. Failure initiation for 12 was by a material compressive strength failure in the comer, 
and an Fj^ = 0.4 at P = 1 .02 P cc is consistent with this observed mode in the sense that delami- 
nation initiation is unlikely. 

The discrepancy between the analysis and experiment for specimen 110 may be due to the 
omission of modeling progressive failures in the laminate. Strain gage data indicated a damage 
event at 18, 750 lbs (P/P cr =1.11, P/P cc = 0.81) for specimen 110, and it is at this load that the 

deviation between analysis and experiment begins to manifest itself in postbuckling on the plot of 
load versus end shortening in Fig. 4. However, first-ply failure criteria did not indicate failure in 
the specimen at 18,750 lbs. Transverse matrix cracks, if present, can influence the neighboring 
interlaminar stress state. For instance, transverse matrix cracks in 90-degree plies can significantly 
reduced the load for delamination initiation at the straight free edge for tensile/compression cou- 
pons (Kim, 1989). Also, Kim (1989) cites difficulties in predicting the onset of delamination 
when it is primarily driven by interlaminar shear stress x^ in combination tensile interlaminar 

normal stress x^. The complexity of the stress fields during postbuckling of the flange, and the 
influence of ply damage on this response, make the prediction of crippling by flange free edge 
delamination a difficult task. 
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Table 1: Experimental data for I-specimens (Bonanni, et al., 1991) 



Specimen 12 Specimen 110 

Average wall thickness t in inches 
Flange width bf in inches 
Web height b w in inches 
Comer radius r in inches 
Specimen length L in inches 

0.0455 0.0833 

1.262 1.263 

1.247 1.25 

0.125 0.125 

7.944 5.916 
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Table 1: Experimental data for I-specimens (Bonanni, et al., 1991) 



Specimen 12 Specimen 110 

Gage length / g in inches 
Buckling load P CT in pounds 
Crippling load P cc in pounds 
Location of damage initiation 

5.944 3.916 

2,786. 16,875. 

10,289. 22,475. 

comer free edge of flange 


Table 2: Material property data used in analysis 


AS4/3502 Graphite/Epoxy 

Specimen 

12 

Specimen 

110 

Fiber direction modulus E\ i in Msi 

17.0 

17.5 

Transverse direction modulus £"22 in Msi 

1.64 

Thickness direction modulus £33 in Msi 

1.64 

Major Poisson’s ratio V 12 

0.3 

Minor Poisson’s ratio V 13 

0.3 

Poisson’s ratio V 23 

0.35 

Shear moduli G \ 2 and G 13 in Msi 

0.87 

Shear modulus G 23 in Msi 

0.49 

0 1 

Longitudinal interlaminar shear strength Z° in ksi 

13.5 

Transverse interlaminar shear strength Z in ksi 

13.5 

Tensile (peel) interlaminar strength Z in ksi 

7.54 

Beam-column material for insert at junctions 



Modulus of elasticity in Msi 

10.0 

13.5 
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Table 3: Maximum value of the dimensionless delamination index F, the corresponding 
interlaminar stresses, and the location of this maximum for specimens 12 and 110. 


ID 

Load 

T7 

Interlaminar stresses for Fj^; psi 

Location along free edge 

P/Pcc a 

r max 

^ZZ 

^zy 

^zx 

x in inches 15 

Interface 0 

12 

1.03 

0.247 

2,674. 

-2,322. 

-4,094. 

4.06 (NP) 

0/90 (3) 


0.92 

0.008 

-893. 

-42. 

839. 

5.48 (BC) 

45/0 (2) 


1.02 

0.400 

-515.5 

1,212. 

8,642. 

6.24 (IP) 

90/0 (4) 

no 

1.03 

0.02 

-1,609. 

1,389 

' 1 

2.80 (BC) 

0/90 (3) 


1.03 

0.73 

-1,961. 

1,438. 

| 

2.03 (IP) 

0/90 (3) 


a. is the crippling load listed in Table L 

b. Letters in parentheses refer to deflection pattern on free edge: NP=>near node point; BC=^near buckle crest; 
IP=>near inflection point. 

c. The number in parentheses corresponds to the interface numbers shown in Fig. 5. 



P 


Fig. 1 The branched-shell finite element model of an l-section stiffener. 
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Fig. 2 The complex model of the l-section stiffener showing a sub-domain containing 
transition and solid elements at the free edge of a flange. 



Fig. 3 Cross section of an l-section stiffener that is fabricated from graphite-epoxy 
unidirectional tape. Unidirectional filler material is required during layup to fill the 
void at the junctions of the flanges and web. 
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End shortening, inches 

Fig. 4 Comparison of the branched shell analysis and experiment for specimen 110 
on the load-end shortening plot. 
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Fig. 5 Distribution of BR20 brick elements through the thickness of a flange, (a) 
Seven elements and six interfaces for specimen 12. (b) Eight elements and seven 
interfaces for specimen 110. 
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Fig. 6 Distribution of the dimensionless delamination index Fin a 0/90 interface for 
specimen 110 at P/P cc =1.03. This is interface 3 in Fig. 5(b). The free edge of the 
flange is at y = 0. 



Fig.7 The distribution of the out-of-plane displacement normalized by the 
flange thickness and the transverse shear resultant Q x along the free edge of the 
flange for specimen 110 from the branched shell model. 
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Fig.8 Distribution of the interlaminar shear stress component tangent to the free edge 
the 0/90 interface of specimen 110 at P * 1.03 P^. 



Fig.9 Distribution of the interlaminar normal stress component in the 0/90 interface o 
specimen 110 at P = 1.03 P cc . 
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II. Pressure Pillowing of an Orthogonally Stiffened Cylindrical Shell 

Background and objectives 

The cabin pressure in a transport aircraft causes about a 10 psi pressure differential across the 
skin. An unstiffened, or monocoque, fuselage would carry this internal pressure load as a shell in 
membrane response, like a pressure vessel. However, internal longitudinal and transverse stiffen- 
ers are necessary to carry maneuver loads, etc. The presence of these internal stiffeners prevents 
the fuselage skin from expanding as a membrane, and the skin bulges, or “pillows”, between the 
stiffeners under the action of the internal pressure. Where the skin is restrained against its expan- 
sion as a membrane along the stiffeners, a bending boundary layer is formed. Also, at the stiffener 
intersection a local concentration of the interacting loads between the stiffeners and the skin 
occurs. 

The objectives of this project are 1) to analyze the concentration of the interacting loads at the 
stiffener intersection, and 2) to study the pillowing of the skin. 

Photocopies of the transparencies used for the presentation at the 34th SDM conference appear in 
the Appendix of this report. These transparencies are the basis of much of the following discus- 
sion, and the page numbers referred to in the text refer to the page numbers labeled on the trans- 
parencies. 

Mathematical model 

An idealized mathematical model is assumed for the semi-monocoque fuselage to study the 
generic characteristics of the response in the vicinity of the stiffener intersection. The model is of 
a very long circular cylindrical shell internally stiffened by identical stringers equally spaced 
around the circumference, and identical frames, or rings, equally spaced along the length as 
shown on page 2. In general, the spacing of the stringers is not the same as the rings. The structure 
is periodic both longitudinally and circumferentially, and the loading is spatially uniform. Conse- 
quendy, a structural repeating unit can be defined whose deformation determines the deformation 
of the entire structure. The repeating unit consists of a portion of the shell wall centered over por- 
tions of the stringer and ring as shown page 3. The radius of the middle surface of the undeformed 
cylinder is denoted by R, stringer spacing by b, and frame spacing is denoted by L. Axial and cir- 
cumferential arc length coordinates are denoted by x and y, respectively. 

The stiffeners are modeled as mathematically one-dimensional structural elements, such that the 
actions transmitted by the stiffeners to the inside surface of the shell wall are represented by dis- 
tributed line loads as shown on page 4. Coordinate £ is an axial coordinate whose origin is cen- 
tered over the stiffener intersection, 0 is the circumferential angle, and £ is the normal coordinate 
to the shell reference surface. It is assumed that the stringer is symmetric about the plane 
through its centroidal axis, and the ring is symmetric about the 0— £ plane through its circular ref- 
erence axis. On the basis of the symmetry about the and 0- axes for the unit, only the interact- 
ing load components tangent and normal to the stiffeners are included in the analysis. The 
interacting load intensities per unit undeformed length along the contact lines are denoted by 
^l(^) for the component tangent to the stringer, for the component normal to the stringer, 

Xy 2 (0) for the component tangent to the ring, and by A-^ 2 (9) for the component normal to the ring. 
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The positive directions for these interacting loads acting on the shell are shown on page 4. The 
purpose of the analysis is to determine these distributed line loads. 

An enclosed volume to contain the pressurized medium is modeled by assuming diaphragms to 
extend from the edges of the repeating unit to the axis of revolution of the cylinder. The dia- 
phragms do not resist the deformation of the repeating unit, but act to transmit loads normal to the 
edges of the repeating unit due to the internal pressure. 

For the linear elastic response of the repeating unit to internal pressure, we use the Ritz method 
and the principle of virtual work augmented by Lagrange multipliers to enforce kinematic con- 
straints between the structural components of the repeating unit. See page 5. The Lagrange multi- 
pliers represent the interacting line loads between the stiffeners and the shell. Displacements are 
separately assumed for the shell, stringer, and ring. The shell is modeled with Sander s theory or 
the Donnell-Mushtari- Vlasov (DMV) quasi-shallow shell theory. 

Displacement and Lagrange multiplier approximations 

The periodic portions of the displacements are represented by truncated Fourier Series, and the 
nonperiodic portions of the displacements due to axial stretching are represented by simple terms 
in For the shell, the displacements of the middle surface are 


M N q£ 

* (i 0) = £ £ “m„ sin ( a m S) cos (P„0) + -jj-. 



m = 1 n = 

0 



M 

N 


v(5,0) 

= i 

£ V m* C0S 

(<x m 4) sin (P„0) , 


m ~ 0 

n = 1 



M 

N 


CD 

= i 

£ w m« C0S 

(a^) cos (p n 0) , 


m = 0 

n = 0 



( 2 ) 

( 3 ) 

( 4 ) 


in which a m = (mti) // and p n = (mz) /O where m and n are non-negative integers, and / = L 
/2. The displacements of the centroidal axis of the stringer are 


M 




M 


«,($> = £ MS) = £ w ,«cos(a m ^) 


m = 0 


( 5 ) 


Parameters qg in the axial displacement of the shell and qj in the axial displacement of the 
stringer represent elongations of each respective element caused by either an axial mechanical 
load or due to closed-end pressure vessel effects. The displacements of the reference circle of the 
ring are 


v r (0) 


N 


£ V r* sin (M) 

n = 1 


M®) 


N 


£ W r« C0S (P« e ) 

n = 0 


( 6 ) 
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( 7 ) 


The distributions of the Lagrange multipliers, or interacting loads, are taken as 


'xl 


M 


X ^i„ sin <°»o 


M 

\i = X \i« cos • 

m = 1 



AT 


Z X >2n Sill ^n e ) 





X V» C0S(P » e) ' 

/I = 0 


( 8 ) 


Note that in the stringer interacting normal load, the term which represents a uniformly dis- 
tributed load, is omitted. Since the stringer as modeled by this analysis is not restrained from rigid 
body motion in the radial direction the resultant normal load on the stringer must vanish; i.e.. 


i 

= 0 , 

-i 


(9) 


and this condition requires Xn o to vanish. 

Discrete equations and their solution 

The discrete displacement vector for the shell is the (3MN+2M+2N +2) x 1 vector 

T 

shell = l®0» ®1» 


( 10 ) 


in which the subvectors are 


fy) “ [<?()’ W 00’ V 01’ w 01> v 0 N’ W 0^ 

I T 


m = 1, M 


(ID 


— W m0’ M ml> V ml’ W ml’ mN • v mN’ 

The (2M + l)x 1 discrete displacement vector for the stringer and the (2N + I) xl vector for the 
ring are 

Kr = tiv tl sl , W gl , ..., U sM , W sM ] T & r = [W^ V rl , W rl> .... V rN , W rN ] T , (12) 


in which the term w s0 for the stringer has been omitted since it does not deform the stringer and its 
conjugate resultant satisfying Eq. (9) is equal to zero. This uniform normal displacement compo- 
nent is determined from the condition that the rigid body displacement of the stringer is the uni- 
form portion of the normal displacement component of the shell at 0 = 0; i.e.. 


N 

X w 0« = w ,o- 

n - 0 

The discrete vectors of the Lagrange multipliers for the stringer and ring, are 


(13) 
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( 14 ) 


^1 = X^ u , •••» ^ 2 “ ^20’ \iv ^C21’ '••\lN’ ^2Jv] ’ 

respectively. 

A schematic of the Lagrange multipliers and the conjugate displacements along the contact 
line between the shell and stringer is shown on page 6. These variables constitute the terms in the 
external virtual work for the shell and stringer. Similar terms are used for the external virtual work 
to represent the contact between the shell and ring. The pressure is treated as a dead load in the 
linear analysis, and contributes to the external virtual work for the shell. The applied axial load P 
is shared between the shell and stringer. To account for load sharing, a Lagrange multiplier Q is 
used to denote that portion of the axial load P earned by the stringer. Thus, P - Q is carried by the 
shell. In the results that are presented below axial load P is prescribed to equal the closed-end 
pressure vessel force that is carried by one repeating unit. 

The discrete system of equations for the linear response of the repeating unit are shown on 
page 7. The first three rows represent equilibrium equations for the shell, stringer, and ring, 
respectively. The fourth row represents the kinematic constraints between the shell and stringer, 
fifth row represents the kinematic constraint between the shell and ring, and sixth row represents 
the fact that the net elongation of the shell and stringer are equal. The right-hand side vector con- 
tains terms determined by the internal pressure. The size of this system for M = N = 24 is shown 
in the table on page 8. 

Numerical example 

Data for the example shown on page 9 are taken from an example solved by Wang and Hsu 
(1985), and the dimensional data is representative of a large transport fuselage. The shell wall 
stiffness matrices given by 



0.664 0.221 

0 ‘ 




262 159 4.33" 

A = 

0.221 0.577 
0 0 

O £1 

d 

ib/in 

5 = 0 

D = 

159 210 4.33 
4.33 4.33 159 


in which A denotes the membrane stiffness matrix, B the coupling matrix, and D denotes the 
bending stiffness matrix of classical lamination theory. 

All the results presented for this example are for an internal pressure p = 10 psi. The Fourier 
Series were truncated at M = N = 24. Since b/R = 2.83°, the shell in this example is shallow and 
the DMV (Donnell-Mushtari- Vlasov) shell theory should be adequate. We found that the numeri- 
cal results using Sanders theory with the rotation about the normal neglected and the numerical 
results using DMV theory were essentially the same. 

Results 

The axial distributions of the shell normal displacement w at three circumferential positions is 
shown on page 1 1 . These displacements for the stiffened shell are all lower than the normal dis- 
placement of 0.229 in. for the membrane shell (no stiffeners). The largest normal displacement of 
the shell occurs midway between the stiffeners, and its minimum occurs at the stiffener intersec- 
tion. The difference in these displacements is a measure of pillowing. Comparison of these results 
to those of Wang and Hsu (1985) (their Fig. 5) show excellent agreement. 
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The axial normal strain distributions on the outside surface and inside surface of the shell are 
shown on page 12. For the data plotted along y = 0, a large bending strain (difference in the sur- 
face strains) gradient is evident as the ring is approached. This axial bending strain gradient at the 
ring midway between the stringers is also a result of panel pillowing. The outside and inside sur- 
face axial strain distributions along the circumference at x = 0 are equal (no bending), and are 
nearly uniform, decreasing slightly as the stringer is approached. Wang and Hsu (1985) show a 
larger decrease in the axial strain as the stringer is approached along the circumferential direction 
midway between the rings (see their Fig. 6). Otherwise, the results shown on page 12 compare 
very well to those presented by Wang and Hsu. 

The interacting normal load intensity acting on the shell due to the contact with the stringer is 
shown on page 13. A positive value means the stringer is forcing the shell radially outward. Since 
the integral of this distribution is zero, the area under this curve is zero. The large spike at the stiff- 
ener intersection is counteracted by small negative mean value over the region outside the inter- 
section. The Fourier Series results for this load suggests a singularity, or non convergence, at the 
stiffener intersection. It was found that for an increasing number of harmonics in the truncated 
series, the peak value at the intersection increased and frequency of the oscillations about the neg- 
ative mean value increased. Wang and Hsu (1985) do not present results for the interacting loads. 

The interacting normal load intensity acting on the shell due to the contact with the ring is 
shown on page 14. A very large negative peak value occurs at the stiffener intersection indicating 
the ring is pulling the shell radially inward. Again this distribution suggests the Fourier Series is 
attempting to model a singularity at the stiffener intersection. 

Peak values of the interacting normal loads acting on the shell at the stiffener intersection are 
plotted versus an increasing number of harmonics (M = N) on page 15. Note that the algebraic 
sum of the two loads is negative, meaning combined effect of both stiffeners is to pull the shell 
radially inward against the internal pressure force. These peak load intensities increase monotoni- 
cally as the number of harmonics increase. Thus, no convergence is apparent in the range of M 
and N shown in the plot. 

The concluding remarks on are given on page 16. The interacting normal load intensities 
exhibit singular behavior at the stiffener intersection, or the stiffener crossing point. Results com- 
pare well to Wang and Hsu’s (1985) results, except for the axial strain decrease in the 0-direction 
at the stringer midway between the rings. 

Investigations of the axial strain distribution in the shell near the stringer 

We tried several methods to investigate if the axial strain decreases in the 9-direction as the 
stringer is approached. First additional terms were added to the displacement expansions, second 
a finite element solution in the axial direction and Fourier Series in the circumferential direction 
was tried, third an analytic method for a cylindrical shell stiffened by stringers only was tried, and 
fourth Wang and Hsu’s (1985) equations were programmed directly. 

Additional displacement terms: The following terms were added to Eqs. (2-4) 

», = ml 1 *™ (W> ”p ■ J*2X C0S (|3 » 9) (1) 

in which q n , / n , and ,y n are additional discrete degrees of freedom. The rationale for these addi- 
tional terms is as follows: A term like £ q\ cos (Pi 9) in the axial displacement permits a circum- 
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ferential variation in the axial strain that is independent of Quadratic terms in \ in the v- 
displacement expansion and the w-displacement expansion are consistent with symmetry about 
the ring, and pemit zero addtional circumferential strain and zero additional in-plane shear strain 
for particular choices of l Q and J n . However, the results with these additional terms showed that 
there were large in-plane shear resultants at % = ± / that did not exist previously. This shear was 
antisymmetrically distributed about the stringer, and zero at 0 = ± 0. The stringer axial load was 
substantially different than the solution without the additional terms shown above. Thus, this 
approach was abandoned. 

Finite element approximation in the x-direction: Starting with the Fourier Series expan- 
sions of the displacements in the circumferential direction, we approximated the fourier series 
coefficients, or generalized coordinates, with finite element interpolations in the axial coordinate. 
We had difficulty in selecting the interpolation functions for the Lagrange multipiers, and could 
not get meaningful results. Thus, we stopped with this approach. 

Analytic solution for stringer-stiffened shell only: We found shell resultants that were inde- 
pendent of axial coordinate q and satisfied equilibrium. Compatibility equations were solved to 
get the displacements. Matching the displacements and force conditions with the stringer resulted 
in a solutions which led us back to spatially uniform resultants and strains. That is, no axial strain 
reduction in the circumferenitial direction near the stringer was obtained. 

Programming of Wang and Hsu’s equations: Finally, we programmed the solution given by 
Wang and Hsu’s paper directly. The result from our coding was that there was no reduction in the 
axial normal strain in the shell as the stringer was approached in the 0-direction. (In Wang and 

Hsu’s paper we calculated coefficients D mn u in Eq. (79) and D n u in Eq. (91) to be zero. Non zero 
values for these coefficients would have given a reduction in the axial strain.) 
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Interacting Loads in an Orthogonally Stiffened Composite 

Cylindrical Shell 
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Long circular cylinder internally stiffened by a regular 
arrangement of frames and stringers 
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Numerical Example 
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present results for them. 
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Axial distributions of the out-of-plane displacement of 
the shell at three circumferential positions 
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Interacting normal load between shell and stringer 
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Interacting normal load between shell and frame 
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